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Abstract 

We present the first rigorous examples of non-singular Hubbard models which ex- 
hibit ferromagnetism at zero temperature. The models are defined in arbitrary di- 
mensions, and are characterized by finite-ranged hoppings, dispersive bands, and finite 
on-site Coulomb interaction U. The picture, which goes back to Heisenberg, that suffi- 
ciently large Coulomb interaction can revert Pauli paramagnetism into ferromagnetism 
has finally been confirmed in concrete examples. 



Introductions: The origin of ferromagnetism has been a mystery in physical science for 
quite a long time [|1] . It was Heisenberg |J who first realized that ferromagnetism is intrinsi- 
cally a quantum many-body effect, and proposed a scenario that spin-independent Coulomb 
interaction and the Pauli exclusion principle generate "exchange interaction" between elec- 
tronic spins. One of the motivations to study the so-called Hubbard model has been to 
establish and understand the generation of ferromagnetism in simplified situations 0, |J. 
Unfortunately, rigorous examples of ferromagnetism (or ferrimagnetism) in the Hubbard 
models have been limited to singular models which have infinitely large Coulomb interaction 
(Nagaoka-Thouless ferromagnetism 0), or in which magnetization is supported by a disper- 
sionless band (Lieb's ferrimagnetism ||, and flat-band ferromagnetism due to Mielke and 
the present author ||). In || IDj , local stability of ferromagnetism in a generic family of 
Hubbard models with nearly-flat bands was proved. 

In the present Letter, we treat a class of Hubbard models in arbitrary dimensions, which 
are non-singular in the sense that they have finite ranged hoppings, dispersive (single- 
electron) bands, and finite Coulomb interaction U. We prove that the models exhibit 
ferromagnetism in their ground states provided that U is sufficiently large. We recall that 
Hubbard models with dispersive bands (like ours) exhibit Pauli paramagnetism when U = 0, 
and remain non-ferromagnetic for sufficiently small U. The appearance of ferromagnetism 
is a purely non-perturbative phenomenon. 

As far as we know, this is the first time that the existence of ferromagnetism is established 
in non-singular itinerant electron systems. We stress that our examples finally provide 
the definite affirmative answer to the long standing fundamental problem; whether spin- 
independent Coulomb interaction can be the origin of ferromagnetism in itinerant electron 
systems ]TT| . See |8], [IT], [T2| for further discussions on ferromagnetism in the Hubbard models. 



1 Phys. Rev. Lett. 75, 4678-4681, 1995. Archived as cond-mat/9509063 
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Figure 1: The one-dimensional lattice and the hopping matrix elements. There are hoppings 
to nearest and next-nearest neighbors, and on-site potential. 

Main results: In order to simplify the discussion, we describe our results in one-dimensional 
models. We discuss models in higher dimensions at the end of the Letter. Let N be an ar- 
bitrary integer, and let A be the set of integers x with \x\ < N. We identify x = —N and 
x = N to regard A as a periodic chain with 2N sites. We denote by S and O the subsets 
of A consisting of even and odd sites (integers), respectively. As usual we denote by c xa , 
c Xi(7 , and n X)CT the creation, the annihilation, and the number operators, respectively, for an 
electron at site x E A with spin cr =|, {. 

We consider the standard Hubbard Hamiltonian 



H 



x,yEA 



i-,a C y 



a + U}^ n x ^n Xji , 
xeA 



where t x>x+1 = t x+1>x = t' for any x E A, t x>x+2 = t x+2jX = t if x E £ , t x>x+2 = t x+2 , x = -s if 
x E O, and t X)X = V if x E O. The remaining elements of t x>y are vanishing. See Figure 1. 
Here s, t, and U are positive parameters fll3| . The parameters t' and V are determined by 
s, t, and another positive parameter A as t' = A(s + 1) and V = (A 2 — 2)(s + 1). Our main 
theorem applies to the case A = v^2, where we have V — 0. We consider the Hilbert space 
with electrons in the system. This corresponds to the quarter-filling of the whole bands, 
or the half-filling of the lower band. 

If we consider the single-electron problem corresponding to the Hamiltonian ([!]), we find 
that the model has two bands with dispersion relations E\(k) = — 2scos2A; — 2(s + 1), and 
£ 2 {k) = 2t cos 2 A; + A 2 (s + 1) with \k\ < tt/2. Note that both the bands have perfect cosine 
dispersions, which is a special feature of the present model There is an energy gap 

A 2 (s + t) between the two bands. 

For a = 1, 2, 3, we define the total spin operators by = J2 x .eA J2a,T=i,i c x a (p^V^)^ c x,t, 
where are the Pauli matrices, and denote the eigenvalues of (S to t) 2 = Sa=i(^tat ) 2 as 
5'tot('S , tot + !)• The maximum possible value of S to t is S max = N/2. 

Let b k a be the creation operator corresponding to the single-electron eigenstate with the 

energy £i(k). Let $ vac be the state with no electrons. The state $f er ro = (Jlk K t)^vac (where 
the lower band is fully filled by up-spin electrons) has the lowest energy among the states 
with Stot = Smax. It is easy to observe that $f er ro is an eigenstate of H with the energy 
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Figure 2: When A = y/2, Theorem I is applicable for t/s and U/s in the shaded region. 
The existence of ferromagnetism is established, for examples, if t > 4.5 s when U — 50 s, or 
t > 2.6 s when U = 100 s. Though the plot was obtained from a numerical calculation in a 
five-site Hubbard model, our theorem guarantees that the model with arbitrary lattice size 
exhibits ferromagnetism. 

Eq = — 2(s + t)N. A simple variational calculation shows that $f err o cannot be a ground 
state of H (and hence the true ground state has Stot < S m3X ) if U < 4s. The main result of 
the present Letter is the following. 

Theorem /—Suppose A > A c = [(2 + V$) 1/2 -2] 1/2 ~ 0.241. If t/s and U/s are sufficiently 
large, the ground states of the Hamiltonian (D) have S tot = Smax, and are non-degenerate 
apart from the (2S tot + l)-fold spin degeneracy. $f err o is one of the ground states. How large 
the parameters should be can be determined by diagonalizing a Hubbard model on a five-site 
chain. (See Figure 2.) 

The present models reduce to the flat-band models studied in H if we set s = 0. Therefore 



we can regard Theorem I as a confirmation in special cases of the previous conjecture [15 



HI that the flat-band ferromagnetism is stable under perturbations. Moreover the strong 



result about the spin-wave excitation proved in || [TO]] also applies to the preset models. 



Theorem II — Suppose that the model parameters satisfy A > A 3 , s/t < p , and K 2 Xt > 
U > v4 3 A 2 s, where A 3 , p , K 2 , and A 3 are positive constants that appear in |10[| . Then the 



spin- wave excitation energy E S y V (k) (i.e., the lowest energy among the states with (N — 1) 
up-spin electrons and one down-spin electron and with crystal momentum k) of the Hubbard 
model ([!]) satisfies 



411 ATI 
F 2 — (sinfc) 2 < E sw (k) -E < F. — ismk) 2 , (2) 



where E is the ground state energy, and F\, F 2 are constants such that F± ~ F 2 ~ 1 if 
A 3> 1, As <S t, and U ^> A 2 s. (See |Tn] for details and a proof.) 



In the parameter region where both Theorems I and II are applicable, we have an ideal 
situation that the global stability of ferromagnetism as well as the appearance of "healthy" 
low-lying excited states are established rigorously. We have rigorously derived a ferromag- 
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netic system with (effective) exchange interaction J ~ 2U/X A , starting from the Hubbard 
models for itinerant electrons! 

It is quite likely that the present models represent (Mott-Hubbard) insulators. We expect 



that the same models with smaller electron numbers describe metallic ferromagnetism [16 



but have no rigorous results in this direction (except for those in the flat-band models 0). 

Proof of Theorem I: For x G A and a =f, |, we define a X)Cr = Xc X)Cr — (—l) x (c x -i^ a + 
c x+i,a), which correspond to the strictly localized basis states used in |8|, pTof] _ The anti- 
commutator {a x CT , a^o-} is A 2 + 2 if x = y, is 1 if \x — y\ — 2, and is vanishing otherwise. By 
using these operators, Hamiltonian ([J) can be written in a compact manner as 

H = {X 2 S - 2t)N - S Y a l,a a x,a + t Y a l,a a x,a + U Y ^,T™*4> ( 3 ) 

xe£ xeo xeA 

f=T)i <T=t,l 

in the sector with N electrons. We further rewrite it as H — (X 2 s — 2t)N + J2xes h x with 
the local Hamiltonian defined as 

h x = s a t,a a x,a + U n xA n xA + Y f I X! <4+r, a<W,<7 + ^ n x+rA n x+rd ] . (4) 

Since [h x , h x+ 2\ ^ 0, it is impossible to diagonalize all h x simultaneously. 

Lemma — Suppose A > A c , and t/s and U/s are sufficiently large. Then the minimum 
eigenvalue of h x (regarded as an operator on the whole Hilbert space) is —(A 2 + 2)s. In 
any of the corresponding eigenstates, there are one, two, or three electrons in the sublattice 
{a; — 2, a; — l, x, x + 1, x + 2}, and these electrons are coupled ferromagnetically. Any eigenstate 
$ with the eigenvalue —(A 2 + 2)s can be written in the form 

$ = 4^ + 4^2, (5) 

with some states $i and $ 2 , and satisfies 

c x+u c x+1 ^ = 0. (6) 

We prove the Lemma in the next part. In what follows, we assume that the model 
parameters satisfy the conditions in the Lemma. The basic strategy of the proof of Theorem 
I is to extend the local ferromagnetism found above into a global ferromagnetism. Special 
characters of the present model makes such an extension possible. 

The Lemma implies h x > -(A 2 + 2)s, and hence H > (X 2 s - 2t)N - N(X 2 + 2)s = 
—2(s + t)N = E . This proves that $f er ro (which has the eigen-energy E ) is a ground state. 

To show the uniqueness of the ground states, we assume $ is a ground state, i.e., = 
E &. Then we have h x <& = —(A 2 + 2)s$ for each x G S, and $ is characterized by the 
Lemma. We note that the collection of states (ILeA °lt) (TIxg-B a li)^vac with arbitrary 
subsets A, B C A such that \A\ + \B\ — N forms a (complete) basis of the A^-electron Hilbert 
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space. Imagine that we expand <£> using this basis. Since (||) holds for any x G £, $ must be 
written in the form 

^ = E^)(II a U.)) $ v-c, (7) 

where a = (a(x)) xe s is a spin configuration with a(x) =|, |, and ip(cr) is a coefficient. Unlike 
in the flat-band models ||, a state of the form (0) is not necessarily an eigenstate of the 
hopping part of H . 

By examining how c x+ i^c x+ \^ acts on (|7|), the condition (§) reduces to 

<y2(cr) = (^(0^+2) for any cr, (8) 

where a x , x +2 is the spin configuration obtained by switching a(x) and a(x + 2) in a. Since 
(|) holds for any x G £, we find that y(cr) = y?(f) whenever Z^e^c^) = J2 X ££ T ( X )- Since 
^ferro is written as ^Wo = const. (flagf a I |)^vac) this means that $ can be written in the 
form $ = EM=o"A/(5'tot) M ^fcrro, where the spin lowering operator is S^ ot = E^eA cl,| c ^.T- 
This proves that $f err o and its SU(2) rotations are the only ground states of H. 

Proof of Lemma: Because of the translation invariance, it suffices to prove the Lemma 
for x — 0. We first diagonalize the hopping part of h (obtained by setting U = 0). We 
express a single-electron state supported on the sublattice A = {—2, —1,0, 1,2} as a five- 
dimensional vector cp = (<£>_2, <P-i, <Pq, <P\, ¥2)- The normalized eigenstates are = (A 2 + 
2)" 1 / 2 (0, -1, A, -1, 0) with the eigenvalue e = -(A 2 +2)s, </? (1) = {2(A 2 +1)}~ 1/2 (A, -1, 0, 1, -A) 
with ei = 0, ^ = {2(A 2 + 2)(A 2 + 3)}- 1 / 2 (-(A 2 + 2), A, 2, A, -(A 2 + 2)) with e 2 = 0, and 
two more with £ 3 = (A 2 + l)t/2 and £4 = (A 2 + 3)t/2. We denote the corresponding creation 
operators as d\ a = J2xeA fx^ c x a- ^ i s crucial to note that a} Q(J = (A 2 + 2) l ^ 2 d\ j(T . 

Since the local Hamiltonian Hq conserves the number of electrons in Ao, we can examine 
its minimum eigenvalue in each sector with a fixed number of electrons in Ao- When there 
are no electrons in A , the only possible eigenvalue of h is > —(A 2 + 2)s. Let f n and e n 
be the minimum eigenvalues of h in the sectors with n electrons in A with the total spin 
(of the n electrons) S^t — n /2, and S^t < n /2, respectively. 

Noting that f n = E^o£m, we find h = f 2 = f 3 = -(A 2 + 2)s and /„ > -(A 2 + 2)s 
for n > 4. Since the corresponding ferromagnetic eigenstates are (Ilmio^t)* (where $ is 
an arbitrary state with no electrons in Ao) or their SU (2) rotations, they are written in the 
desired form (|5]), and satisfy (|6|). Therefore, in order to prove the Lemma, it suffices to show 

e n > -(A 2 + 2)s for any n = 2, 3, . . . , 8. (9) 

Since the condition (|S]) only involves eigenvalues of a finite system, it can be checked by 
numerically diagonalizing finite dimensional matrices for given values of A, s, t and U. We 
can thus construct a computer aided proof that our Hubbard model exhibits ferromagnetism. 
Figure 2 summarizes the result of a preliminary analysis in this direction. 

Let us prove (H) in a range of parameters without using computers. Let e s 2 yra (resp., ejf) 
be the minimum eigenvalue of ho in the sector with two electrons in Aq forming spin-singlet 
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states which is symmetric (resp., antisymmetric) under the spatial reflection x — > —x. Let 
us evaluate e^ 7 ™. In the limit t | oo, a spin-singlet state with two electrons in the symmetric 
sector which has finite expectation value of ho is written as 

* = |«4,t<! + -^(4A,i + 4 T 4,i) +7 4,t4i + 5 4A} ^ ( 10 ) 

where $ is any state with no electrons in Ao- The expectation value of the hopping part of 
ho in this state is given by (h^) = /i£ op *)/(*, *) = A[-2(A 2 + 2)s|a| 2 - (A 2 + 2)s\(3\% 
where A = (\a\ 2 + \fl\ 2 + |7| 2 + |5| 2 ) -1 . If we further let U \ oo, a finite energy state must 
also satisfy Cq^Cq^I = and ci^ci^ = 0. These conditions lead us to the constraints 



A 2 + 2 (A2 + 2 )v / FT3^ (A 2 + 2)(A 2 + 3) 

' r«~ 7^ ■ ^ ■ o fl+ oM, ■ £vi ■ o^ 7 + T77T7 L ^^ = 0• (") 



A 2 + 2 (A2 + 2)v / A T T3^ 2(A 2 + 2)(A 2 + 3) ' 2(A 2 + 1) 

We denote the desired minimum eigenvalue e^ m for t = U = oo as e^™. To get e^™, we 
minimize the energy expectation value (/io° p ) with respect to the constraints (|Tl~l) . The rest 
is a tedious but straightforward estimate. Eliminating (3 from (|llD , we get a + 7/ (A 2 + 3) + 
S/(\ 2 + 1) = 0, which implies | 7 | 2 + \S\ 2 > f(\)\a\ 2 with /(A) = {(A 2 + l)" 2 + (A 2 + 3)- 2 }- 1 . 
By substituting this bound into (ho) + (A 2 + 2)s = A(X 2 + 2)s(-|a| 2 + | -y | 2 + \5\ 2 ), we get 
+ (A 2 + 2)s > A(\ 2 + 2)s|a| 2 {/(A) - 1}. Noting that a ^ in the minimizer, the 
condition /(A) > 1 (which is equivalent to A > A c ) implies > —(A 2 + 2)s. Since e^ m is 
a continuous function of t and U, this proves that e^ m > —(A 2 + 2)s for sufficiently large t 
and U. 

By repeating the similar (but easier) variational analysis, we find that eff = — (A 2 /3)s > 
— (A 2 + 2)s, and e n = 00 for n > 3 when t = U = 00. This implies that the desired condition 
(0) holds for A > A c and sufficiently large t and U. The Lemma has been proved. 



Models in higher dimensions: Models in higher dimensions can be constructed and 



analyzed in quite the same spirit |12l . Take, for example, the flat-band models studied in [f 
(V and M correspond to £ and O of the present paper, respectively.) For x G V we let a a 



as m 

as 



For x = m(v, w) G M, we let a a 



XC X „ ~\~ C v a ~\~ 



We define the Hamiltonian 



H = -s a l,a a x,a + * E ' 

x&V x£M 

t=T4 <r=T,J. 



;<r + U ^ ^.T^ar.l, 



(12) 



which again contains next-nearest neighbor hoppings and on-site potentials. should be 
compared with @. By a straightforward extension of the present method, we can prove 
that the ground states of the model with \V\ -electrons exhibit ferromagnetism when A, t/s, 



and UJ s are sufficiently large [Tj| . 
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Note added (August 1997): The proof of Theorem I has been considerably improved. 
Now the condition A > A c has been replaced simply by A 7^ 0. Full details will appear in 



12], which is still under preparation for the moment. 
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